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ABSTRACT 


•-^The  ostensible  purpose  of  a mathematical  progr^UMning  model  is  to 
optimize  a stipulated  objective  function  subject  to  stipulated  constraints. 
But  its  true  purpose,  at  least  in  strategic  applications  as  every  ex- 
perienced practitioner  should  knew,  is  to  help  develop  insights  into 
system  behavior  which  in  turn  can  be  used  to  guide  the  development  of 
effective  plans  and  decisions.  Such  insights  are  seldom  evident  from  the 
output  of  an  optimization  run.  One  must  know  not  only  whatthe  optimal 


solution  is  for  a given  set  of  inpat  data,  but  also  why.^^  ~^e  desired  insights 
usually  have  more  to  do  with  the  "why"  than  the  "whaf'.SiThis  paper  advocates 
the  use  of  highly  simplified  analytic  models  to  help  explain  the  ■'^rtiys^be- 
hind  the  solutions  of  conventional  mathematical  programming  model^^  A 
metiiodological  approach  is  described  which  permits  the  development  of 
richer  insights  than  would  otherwise  be  possibl^^^ This  approach  is  illus- 
trated with  reference  to  a facility  location  study  carried  out  recently  for 
a consumer  products  manufacturer. 
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I.  A METHODOLOGICAL  PROPOSAL 

Mathematical  programming  models  are  fine  as  far  as  they  go.  The 
trouble  is  that  they  seldom  go  far  enough  due  to  a serious  inherent  limi- 
tation: they  can  deliver  solution  for  a given  set  of  input  data, 

but  they  do  not  explain  WH\  . , -ion  is  what  it  is.  Numerical  models 
involving  an  iterative  solution  procedure  engender  a 'humbers  in-nuitibers 
out"  user  attitude  that  tends  to  inhibit  the  development  of  insights  due  to 
a false  sense  of  (numerical)  security. 

Yet  it  is  very  in^jortant  to  develop  fundamental  insights  into  the 
reasons  why  an  optimal  solution  is  what  it  is,  particularly  for  strategic 
applications  of  mathematical  programming.  One  reason  is  that  few  if  any 
applications  lead  to  a single  perfect  numerical  model  whose  solution  is 
directly  tremsla table  into  practical  action.  Rather,  there  is  a family  of 
imperfect  numerical  models  reflecting  alternative  assumptions,  objectives, 
and  data  estimates;  an  understanding  of  the  solution  behavior  of  the  whole 
family  is  needed  in  order  t o fully  support  the  development  of  an  appro- 
priate plan  of  action. 

Another  reason  for  the  importance  of  iasights  into  the  determinants  of 
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an  optimal  solution  is  that  they  help  to  overcome  the  serious  validation/ 
credibility  obstacles  so  often  present  in  practical  applications.  How  can 
one  be  convinced  a model  is  a useful  representation  of  the  real  system? 

And  how  can  the  ond-usor  of  a model  — usually  a managerial  or  political 
figure  rather  than  a technical  person  — be  percuaded  to  use  the  model  as 
a problem-solving  aid?  The  answer  to  both  questions  is,  in  substantial 
part,  that  purely  numerical  results  must  be  supplemented  by  intuitively 
reasonable  explanations  as  to  why  these  results  are  as  they  are. 

Otherwise  the  validity  of  a model  can  only  be  taken  as  an  act  of  faith 
a:id  the  end-user  will  be  inclined  to  revert  tc  intuition  or  some  other 
more  secure  mode  of  analysis. 

Where  will  the  desired  insights  come  from  to  illuminate  numerical  mathe- 
matical programming  models?  The  approach  advocated  here  is  that  they  should 
come  from  simplified  auxiliary  models  that  are  both  intuitively  plausible  and 
solvable  in  closed  form  or  by  simple  arithmetic.  The  solution 
behavior  of  a well-chosen  auxiliary  model  should  be  vastly  more  transparent 
tlian  that  of  the  full  mathematical  programming  model,  yet  it  should  yield 
fairly  good  predictions  of  the  general  solution  characteristics  of  the  full 
model. 

A general  methodological  approach  is  as  follows: 

1.  Reduce  the  level  of  detail  and  complexity  of  the  full  mathematical 
programming  model  until  it  can  be  solved  in  closed  form  or  by 
simple  arithmetic.  Call  this  an  aaxiJUMuj  modeJL. 


2.  Derive  from  the  auxiliary  model  a set  of  tentative  hypotheses 


concerning  the  general  behavior  of  the  solution  of  the  full 
Rodel  — the  cost  tradeoffs  determining  the  optimal  solution 
for  a given  set  of  data,  t)ie  nature  of  the  induced  change  in 
the  optimal  solution  as  certain  data  are  changed  parametrically, 
and  so  on. 

3.  Generate  specific  predictions  from  the  tentative  hypotheses  and 
test  these  numerically  using  the  full  model. 

4.  To  the  extent  that  the  numerical  tests  confirm  (i.e.,  do  not  contra- 
dict) the  tentative  hypotheses,  tedee  these  hypotheses  as  a 
conceptual  framework  for  understcuiding  cmd  interpreting  the 
numerical  results  provided  by  the  full  model. 

The  use  of  relatively  simple  models  in  conjunction  with  conplex  models 
is  not  a new  idea.  "Pilot"  or  otherwise  scaled-down  versions  of  full-scale 
models  have  long  been  used  as  a practical  evolutionary  approach  to  model- 
building. Since  solving  them  still  requires  elaborate  algorithmic  machinery, 
however,  they  are  generally  unsuited  to  the  purposes  we  have  outlined  for 
auxiliary  models.  So-called  "repro-models"  [8]  have  been  proposed  as  simpli- 
fied input/output  approximations  to  conplex  models.  They  are  mathematical 
constructs  which  lack  the  ejq)lanatory  properties  essential  for  an  auxiliary 
model.  Much  closer  to  the  auxiliary  model  concept  is  the  use  of  sinple 
cuialytical  models  that  are  validated  and  perhaps  calibrated  with  the  help  of 
more  complex  numerical  models.  Ignall,  Kolesar  emd  Walker  give  several 
splendid  exeunples  of  this  approach  in  the  context  of  emergency  service  de- 
ployment [5]  . It  so  happens  that  their  full-scale  models  are  of  the  simu- 
lation rather  them  mathematical  programming  variety.  The  main  difference 
between  this  and  our  approach  to  auxiliary  models  is  apparently  one  of  re- 
versed objectives:  they  want  full-scale  models  to  support  the  use  of  simple 

models,  whereas  we  want  simple  models  to  support  the  use  of  full-scale 
models.  This  reversal  of  objectives  is  probably  a consequence  of  the  fact 
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that  they  were  primarily  concerned  with  tactical  operational  applications 
where  elaborate  simulations  may  be  cumbersome,  whereas  we  are  primarily  con- 
cerned with  strategic  planning  applications  where  the  resources  needed  to 
support  elciborate  mathematical  progr2uraning  models  may  be  more  readily 
available.  In  any  case,  I definitely  intend  auxiliary  models  to  supplement 


full-scale  models,  not  to  replace  them. 

The  next  section  summarizes  my  experience  with  auxiliary  models  in  the 
context  of  a facility  location  study  carried  out  recently  for  a consumer  pro 
ducts  manufacturer.  Three  different  auxiliary  models  were  developed  in  am 
effort  to  understand  and  eoqilaln  the  solution  behavior  of  the  large  integer 
programming  model  (over  1400  0-1  variables  amd  250  constraints). 


II.  APPLICATION  IN  THE  CaiTEXT  OF  FACILITY  LOCATION 


A consumer  products  manufacturer  has  a single  plant  in  California 
making  a range  of  products  which,  for  the  purposes  of  a warehouse  location 
study,  can  be  treated  as  a single  product  group.  Goods  are  distributed 
nationally  through  full-line  warehouses  to  customers  grouped  into  about 
130  customer  zones.  There  are  about  65  permissible  warehouse  locations, 
tlie  best  subset  of  which  is  to  be  selected.  Business  policy  requires  that 
each  customer  zone  must  be  single-sourced,  that  is,  it  may  not  receive 
goods  from  more  than  one  warehouse.  The  main  objective  is  to  decide 
which  warehouses  to  use  and  to  design  their  service  areas  so  as  to  mini- 
mize the  sum  of  all  relevant  costs  (freight  in,  fixed  and  variable  ware- 
house-related costs,  and  freight  out)  subject  to  the  constraints  that 
forecast  annual  demands  will  be  met,  customer  service  standards  will  be 
maintained,  and  warehouse  aggregate  throughput  limits  (both  lower  and 
upper)  will  be  honored. 

The  problem  as  just  described  is  a familiar  one  in  the  literature  on 
facility  location  (e.g.,  [3]).  It  has  been  simplified  for  the  purposes  of 
this  illustration  by  omitting  mention  of  several  complexities  treated  in 


the  actual  study,  including:  plemt  direct  shipments,  economies  of  scale 

on  inbound  freight  ctnd  warehouse  throughput  costs,  the  conplexities  of 
customer  will-calls  and  delivery  consolidation  services,  and  the  need  to 
choose  between  alternative  inventory  stockage  policies. 

I shall  describe  in  turn  three  of  the  auxiliary  models  that  have  been 
used  in  conjunction  with  the  big  integer  progreunning  model. 
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A.  The  First  Auxiliary  Model 

The  first  auxiliary  model  ia  based  on  one  of  the  simplest  related  pro- 
blems to  be  found  in  the  location  economics  literature.  It  meikes  the  follow- 
ing simplifying  assumptions: 


Al.  Demand  is  uniformly  distributed  on  the  plane  with  a density  of 
2 

P (CWX/mi  ) . 

A2.  All  warehouses  are  identical  and  arbitrarily  relocatable. 

A3.  The  supply  cost  for  each  warehouse  is  s ($/CWT)  regardless  of 


its  location. 


A4.  The  fixed  cost  of  each  warehouse  is  f ($). 


A5.  The  variable  throughput  cost  of  each  warehouse  is  v ($/CWT)  . 

A6.  The  o itbound  freight  rate  for  each  warehouse  is  t ($/CWT-mi.). 

A7.  There  are  no  throughput  limits  for  the  warehouses. 

How  many  warehouses  should  there  be,  where,  and  with  what  service 

2 . . 

areas  to  cover  a total  area  A (mi  ) of  unspecified  shape  at  minimum  total 
cost?  The  answer  (e.g.,  [1],[6])  is  as  follows:  there  should  be 


n*  = A (P  t/f) 
3.05 


warehouses—'  , each  at  the  center  of  a circular  service  area  comprising 


A*  = 3.05  (P  t/f) 


square  miles,  located  anywhere  on  the  plane  so  long  as  their  service 


areas  do  not  overlap. 

V Actually,  n*  as  given  can  be  fractional;  in  this  case  it  is  necessary 
to  examine  the  total  cost  expression  for  the  two  integers  on  either 
side  of  n*.  We  shall  neglect  tliis  detail  here  because  n*  is  quite 
large  in  this  particular  application. 


* 


It  is  also  known  that  this  optimal  solution  is  remarkably  insensitive 


to  the  exact  geometry  of  tlie  warehouse  service  areas.  Hexagonal  rather 
than  circular  service  areas  cliange  n*  by  only  about  + 0.2%,  and  even 
square  service  areas  would  change  n*  by  only  about  + 1.1%.  This  insensi- 
tivity, along  with  the  arbitrariness  of  location  (one  need  only  avoid  over- 
lapping service  areas),  enables  one  to  position  the  warehouses  euid  adjust 
their  service  area  shapes  so  as  to  give  a fair  approximation  to  the  shape 
of  the  United  States.  The  degree  of  suboptimality  thereby  incurred  is 
likely  to  be  very  small  if  n*  is  more  than  a dozen  or  so  , in  which  case 
(1)  and  (2)  give  a very  nearly  optimal  solution  for  the  more  realistic 
situation  where  the  demand  to  be  covered  is  uniform  over  the  continental 
U.S. 


One  cam  readily  answer  almost  any  "why"  question  associated  with  this 
simple  model.  The  answer  at  a mathematical  level  is  available  if 
one  understands  the  mathematical  derivations  supporting  (1)  auid  (2) . 

This  is  in  stark  contrast  to  the  clumsiness  of  a large  numerical  integer 
programming  model.  Effective  communication  of  these  insights  at  a 
managerial  level,  however,  requires  that  a clear  rationale  be  given  for 
(1)  and  (2)  in  terms  fcimiliar  to  management.  The  key  is  to  interpret  A* 
as  that  circular  service  area  size,  with  the  warehouse  at  the  center, 
which  minimizes  the  average  unit  cost  of  satisfying  all  demand  within 

the  area.  Supply  and  variable  warehouse  throughput  costs  can  be 
ignored  because  their  cost  contributions  are  independent  of  warehouse  loca- 
tion and  service  area  shape.  That  leaves  fixed  costs,  which  clearly  are  in 
versely  proportional  on  a unit  basis  to  service  area  size,  and  outbound 
freight  costs,  which  increase  proportionally  on  a unit  basis  with  the 
square  root  of  service  area  size  (this  follows  from  purely  dimensional 
considerations).  For  a service  area  size  less  (more)  than  A*  as  given  by 


Fixed  and  Outbound  Freight  Costs  Per  square  mile 
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(2),  die  unit  fixed  costs  are  decreasing  faster  (slower)  than  the  unit  out- 
bound freight  costs  are  increasing.  A graph  li)te  Figure  1 brings  these 
words  alive  for  manager  and  management  scientist  alike.  Here  we  have  used 
average  values  of  p,  f,  t and  A obtained  from  the  detailed  data  prepared 
originally  to  support  the  integer  programming  model  for  the  consumer  pro- 
ducts manufacturer. 


If  a warehouse  with  a circular  service  area  of  size  A*  gives  the  minimum 
cost  per  square  mile  of  meeting  demand,  then  A/A*  surely  gives  the  number  of 
warehouses  needed  to  meet  demand  over  a total  area  A in  the  most  economical 
way  (presuming  A/A*  is  an  integer).  Indeed,  n*  in  (1)  can  be  rewritten 
equivalently  as  n*  = A/A*. 

Such  an  explanation  gives  a quite  satisfactory  intuitive  understanding 
for  (1)  and  (2).  The  tradeoff  between  fixed  costs  and  outbound  freight 
costs  is  clearly  revealed  as  the  factor  which  controls  the  optimal  solution. 
Hopefully  this  insight,  even  though  gained  for  but  a grossly  simplified 
version  of  the  problem,  will  contribute  to  an  understanding  of  the  behavior 
of  the  real  system.  One  way  to  test  the  validity  of  this  reasoning  is 
to  generate  sensitivity  analysis  predictions  using  the  auxiliary  model  and 
then  to  test  these  using  the  integer  programming  model  (refer  to  Steps  2-4 
of  the  general  methodological  approach  given  at  the  end  of  Sec.  I) . This  will 
now  be  done. 

Relation  (1)  implies  an  easily  calculated  percentage  change  in  n*  for 

any  given  percentage  change  in  p , t or  f.  For  instance,  increasing  p or  t by 

2/3 

5%  implies  that  n*  will  increase  by  3.3»  (since  1.05  » 1.033).  These 
sensitivity  relationships  were  examined  empirically  using  the  integer  pro- 
gramming code  for  several  sets  of  data.  The  results,  given  in  Table  1,  in- 
dicate that  the  sensitivity  predictions  are  svyrisingly  accurate  for  small 
changes  in  p , t amd  f . 

The  successes  shown  in  Table  1 add  to  our  confidence  in  the  auxiliary 
model  as  a source  of  insight  into  the  behavior  of  the  real  system.  They 
suggest  that  the  explanations  for  the  auxiliary  model's  sensitivity  be- 
havior, which  can  be  spelled  out  clearly  in  full  detail,  may  also  be 
operative  for  the  integer  programming  model  and  hence  (hopefully)  for  the 
real  system. 


Integer  Proqranuning  Model 


Optimal  n 
before  change 


Predicted  n 
after  change^ 


24 

24 

24.8 

24 

25 

25.1 

26 

28 

28.2 

24 

25 

24.9 

Table  1 

Increase  p by  5% 
Increase  t by  7% 
Increase  t by  13% 
Decrease  f by  5% 


Sensitivity  Analysis  with  the  First  Auxiliary  Model 
(data  for  consumer  products  memufacturer) 

Another  use  to  which  the  auxiliary  model  can  be  put  is  to  examine  the 
conventional  wisdom  quoted  so  often  by  practitioners  of  the  art  of  facility 
location  analysis:  that  total  system  costs  are  quite  insensitive  to  moderate 

departures'  from  the  optimal  number  of  warehouses.  The  usual  way  in  which 
this  is  demonstrated  is  to  calculate  the  total  costs  associated  with  several 
manually  configured  systems,  each  with  a different  number  of  warehouses  and 
with  service  areas  drawn  as  well  as  manual  methods  will  allow  (e.g.,  Magee 
[7]).  The  shortcoming  of  this  approach  is  that  it  does  not  use  a least  cost 
system  associated  with  a given  number  of  warehouses;  the  warehouse  locations 
and  their  service  areas  are  not  necessarily  optimal.  This  shortcoming  cein  be 
overcome  in  the  context  of  the  auxiliary  model.  Given  that  there  will  be  n 
warehouses,  the  optimal  system  under  the  assumptions  of  the  auxiliary  model 
has  total  cost 


(3)  TC*(n)  = pA(v  +s)  + nf  + .377  ptA5'n"2 

The  first  term  is  independent  of  n and  can  be  ignored.  The  second  (fixed 
cost)  and  third  (outbound  freight  cost)  terms  emd  their  sum  are  graphed 
in  Figure  2 using  the  S2une  average  values  for  p,  A,  f and  t as  were  used 

2/  The  predictions  were  computed  as: 


n*  from  (1)  after  data  chemges 
n*  from  (1)  before  data  chcinge 


n*  from  IP  model  1 
before  data  change. 


1 


for  Figure  1.  The  ordinate  is  given  in  relative  costs,  with  1.0  being  the 
fixed  plus  outbound  freight  costs  associated  with  n*.  It  is  evident  from 


Figure  7.  t)iat  there  is  indeed  considereible  latitude  to  depart  from  n* 
without  Incurring  a substantial  cost  penalty.  For  instance,  n can  range 
between  about  33  and  51  and  yet  keep  fixed  plus  outbound  freight  costs 
within  1%  of  optimum.  The  explanation  for  this  relative  insensitivity  is 
evident  from  the  graphs  of  the  component  cost  functions  in  Figure  2. 


fixed  plus 


Number  of  Warehouses 
Figure  2 


Influence  of  the  Number  of  Warehouses  on  an  Otherwise  Optimally 
Configured  System  According  to  the  First  Auxiliary  Model 
(data  for  consumer  products  m2mufacturer) 
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An  obvious  limitation  of  this  first  auxiliary  model  is  that  assump- 
tion A3  is  overly  stringent.  Indeed,  A3  virtually  ignores  the  economic 
impact  of  plant  location  and  inbound  transportation.  The  next  auxiliary 
model  is  designed  to  overcome  this  limitation. 


B.  The  Second  Auxiliary  Model 

In  order  to  incorporate  the  influence  of  inbound  transportation  eco- 
n mics  it  seems  necessary  to  replace  assumption  A3  by  a version  in 
which  inbound  freight  costs  are  given  in  terms  of  a $/cwT-mi  rate  instead 


of  a flat  $/CWT.  If  all  other  assumptions  remain  the  same,  however,  this 
change  introduces  an  exceedingly  messy  integral.  I have  therefore  taken 
the  more  expedient  course  of  working  out  a one-dimensional  version  of  the 
first  auxiliary  model.  The  assunptions  are: 


Ala.  Demand  is  uniformly  distributed  on  a straight  line  of  length  L 
(mi.)  with  a density  p (CWT/mi). 

A3a.  There  is  a single  plant  at  one  end  of  the  line,  supplying  each 
warehouse  at  cin  inbound  freight  rate  r ($/OTT-mi) , where 


A2,  A4-A7  as  before. 


How  many  warehouses  should  there  be,  where,  and  with  what  service 
areas  to  cover  the  demand  over  the  entire  line  at  minimum  total  cost?  It 


can  be  shown  that  there  should  be 


n*  = f b ( P t/f) 
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wartihouses-^,  each  with  a service  area  of  size  i * L/n*  miles,  configured 


as  shown  in  Figure  2.  Notice  that  each  warehouse  is  displaced  from  the 
center  of  its  service  area  in  ttie  direction  of  the  plant  by  an  amount  pro- 
portional to  r/t.  The  reader  may  find  it  interesting  to  note  the  simi- 
larities and  dissimilarities  between  e)q>ressions  (1)  and  (4)  . 


Figure  3 

Optimal  Service  Area  Configuration  for  the  Second  Auxiliary  Model 


A rationale  to  eo^lain  in  simple  managerial  terms  why  the  above  con- 
stitutes an  optimal  solution  ie  not  given  here,  but  can  be  worked  out  along 
the  lines  of  the  rationale  r^’esented  for  the  first  auxiliary  model. 


The  reader  may  be  curious  as  to  how  close  n*  from  (4)  comes  to  the 
value  obtained  from  the  integer  programming  model:  n*  from  (4)  is  26, 


w!  lie  the  integer  programming  optimum  was  24.  Surprisingly,  this  is  better 


than  the  estimate  given  by  the  first  auxiliary  model  (for  which  n*  was  Al—  ) 
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This  ntodel  has  some  remarkable  in^lications  concerning  the  in- 
fluence of  plant  location  and  inbound  transportation  costs  on  the  optimal 
warehouse  configuration.  We  see  from  (4)  emd  Figure  3 that  this  influence 
enters  only  via  the  ratio  r/t.  This  ratio  was  .082  for  the  consumer  pro- 
ducts manufacturer,  a reflection  of  the  fact  r was  based  mostly  on  rail 
carload  rates  but  with  some  full  truckload  shipments,  while  t was  based 
mostly  on  much  higher  less-than-truckload  rates.  The  consequences  of  such  a 
low  r/t  ratio  are  surprising.  According  to  (4) , n*  is  essentially  the  same 


as  though  r were  0,  since  /l  - (.082)  = .9966  = 1.  And  each  warehouse 

is  displaced  from  the  center  of  its  service  area  toward  the  plcint  in 
California  by  a mere  8.2%  of  the  half-width  of  each  service  area  (less 
than  5 miles) . 

In  other  words,  this  model  indicates  that  for  firms  with  such  a low 
r/t  ratio  the  influence  of  inbound  transportation  economics  and  plant 
location  on  the  optimal  warehouse  configuration  is  extremely  small.  In 
such  cases  one  may  as  well  neglect  inbound  transportation  altogether,  as 
the  first  auxiliary  model  essentially  did. 


An  attempt  was  made  to  confirm  empirically  the  predicted  insensitivity 
with  respect  to  r.  According  to  (4),  the  r/t  ratio  could  vary  from  0 to 


more  than  0.2  without  significantly  altering  n*  . In  fact,  however, 
scaling  all  inbound  freight  rates  by  a factor  of  j (which  reduces  r/t  to 
.041)  causes  the  number  of  warehouses  in  the  optimal  solution  of  the  in- 
teger programming  model  to  increase  from  24  to  28.  This  observation  was  in 
clear  disagreement  with  the  model's  prediction.  It  seemed  at  first  to  cast 
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serious  doubt  on  the  usefulness  of  the  second  auxiliary  model,  but  further 
reflection  revealed  tliat  the  empirical  test  could  have  been  prejudiced  by 
a confounding  effect:  factoring  all  inbound  freight  rates  by  1/2  not  only 

halved  the  average  as  desired,  but  also  halved  the  deviation  about  the 
average  for  individual  rates.  The  latter  effect,  which  cemnot  be  reflected 
in  the  second  auxiliary  model,  can  be  eliminated  simply  by  reducing  the 
inbound  freight  rates  in  a subtractive  rather  than  multiplicative  fashion. 
More  precisely,  subtract  half  the  average  inbound  rate  from  each  individual 
rate  instead  of  factoring  each  individual  rate  by  1/2.  This  halves  the 
average  rate  as  before  but  leaves  unchanged  the  deviation  of  individual 
rates  about  the  average.  The  integer  programming  model  was  rerun  using  the 
subtractive  reductions.  As  hoped,  this  integer  programming  solution  ex- 
hibited the  relative  insensitivity  predicted  by  the  second  auxiliary 
model:  the  optimal  number  of  warehouses  increased  by  just  one. 

Thus  we  are  led  by  the  second  auxiliary  model  to  an  importamt  in- 
sight concerning  the  behavior  of  the  full  integer  programming  model:  ihz 

optunaZ  iolwUon  Zi  quJXz  Zn&zn&itivz  to  dzviation-pKziZAving  chcmgzi  in 
thz  mzan  inbound  ^fieight  Katz  whin  thi  K/t  Kotio  it  (Technically,  of 

course,  this  "insight"  should  more  properly  be  called  a tentative  hypothesis 
for  which  there  is  some  analytical  and  empirical  evidence.) 

The  influence  of  deviation-altering  chaunges  in  the  Inbound  rates,  as 
well  as  other  warehouse-specific  data  sensitivities,  can  be  studied  via  a 
third  auxiliary  model. 


i 

f 

r 
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C.  The  Third  Auxiliary  Model 

This  auxiliary  model  Is  Identical  with  the  first  except  that  It  drops 
the  assumption  that  costs  are  identical  for  all  warehouses.  Let  k index 
the  candidate  warehouses. 

Al,  A7  as  before. 

A2a.  All  warehouses  are  arbitrarily  relocatable. 

A3b.  The  supply  cost  for  warehouse  k is  Sj^($/CWT)  regardless  of 
its  location. 

A4a.  The  fixed  cost  of  warehouse  k is  f,  ($)  . 

k 

A5a.  The  variable  throughput  cost  of  warehouse  k is  Vj^($/CWT). 

A6a.  The  outbound  freight  rate  for  warehouse  k is  tj^ ( $/CWT-mi ) . 

It  is  not  possible  in  general  to  derive  a closed  form  optimal  solu- 
tion along  the  lines  of  (l)-(4).  There  is,  however,  a simple  non- 
iterative procedure  for  calculating  how  many  warehouses  there  should  be  emd 
with  what  size  service  areas  (emd  even  which  ones  should  be  selected) . It 
will  be  convenient  to  work  with  the  demand  Dj^(CWT)  covered  by  the  circular 
service  area  centered  on  warehouse  k,  rather  than  in  terms  of  the  actual 

area  A,  covered  (D,  = p A,  ) , Instead  of  having  to  cover  a total  area  A we 
k k k 

shall,  equivalently,  cover  a total  demand  D (D  = p A) . 

The  procedure  for  solving  the  problem  is  as  follows.  Its  validity 
should  be  evident  from  a rationale  similar  to  the  one  given  previously 


.'J 


j 


for  the  first  auxiliary  model. 
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Step  1 Calculate  for  each  k the  value  D*  which  minimizes 

_ 1 

“k  “k 

where  ATCj^(Dj^)  measures  the  average  total  cost  in  $/ClVT  associated 

with  warehouse  k when  it  has  a circular  service  area  covering  an 
amount  Dj_  of  demand.  Denote  ATCj^(D*)  by  AfC*  . 

Step  2 Sort  on  ATC*  and  reii.dox  the  warehouses  so  that 
(6)  ATC*  < ATC*  < ATC*  < ... 

K 

Step  3 Plot  the  partial  sums  E D*.  ATC,*  against  the  partial  sums 

k=l 

K 

E D*  for  K = 1,2,...  emd  mark  D on  the  abcissa.  See  Figure  4. 
k=l 

step  4 If  D corresponds  with  one  of  the  partial  sums  of  D*,  say  the 

K then  K is  the  optimal  number  of  warehouses  and  the  particular 

o o 

ones  which  should  be  used  are  indexed  1 <k  <K  . The  k^'  should  have 

- - o 

a circular  service  area  covering  D*  CWT  of  demand,  and  they  can  he 
located  anywhere  in  the  plane  so  long  as  their  service  areas  do  not 
overlap.  In  Figure  4,  = 4. 

If  D does  not  correspond  exactly  to  one  of  the  partial  sums  of 
D*|^,  then  cui  optimal  solution  is  no  longer  immediately  obvious.  However, 
unless  D is  relatively  small,  it  is  adequate  for  our  purposes  to  allow  K 

o 

to  be  fractional,  just  as  n*  can  be  fractional  in  (1)  and  (4).  Another 

essential  comment  on  Step  4 is  that,  as  expledned  in  the 

context  of  the  first  auxiliary  model,  only  a small  degree  of  suboptimality 
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is  incurred  if  K is  larger  them  a dozen  or  so  and  the  warehouses  are 
o 


placed  as  necessary  with  reshaped  service  eureas  so  as  to  cover  the 
continental  U.S. 


< 


♦ a: 
n 


II 
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Figure  4 

Hypothetical  Illustration  of  the  Solution  Procedure 
for  the  Third  Auxiliary  Model 


This  model  is  so  closely  related  to  the  first  auxiliary  model  that  a 
very  similar  rationale  in  mautagerial  terms  can  be  given  for  the  "why"  behind 
its  solution.  For  Instemce,  the  cost  relationships  shown  in  Figure  1 still 
apply. 

How  well  does  this  model  predict  the  behavior  of  the  integer  prograirming 
model  for  the  consumer  products  manufacturer?  Table  2 gives  conparative  re- 
sults in  terms  of  n*  and  in  terms  of  a cost  breeikdown.  The  auxiliary  model 
was  applied  twice,  once  with  the  procedure  as  stated  and  once  with  a simple 
manual  modification  which  took  account  of  the  fact  that  closely  adjacent  ware- 
houses cannot  be  open  simultweously  if  their  service  areas  would  obviously  over 
lap.  The  agreement  in  both  cases  is  fairly  good.  Moreover,  the  auxiliary  model 


Optirtal  Number  of  Warehouses 
Optimal  Costs(percentagewis8) 
Fixed 

Inbound  + Var.  Thruput 
Outbound  Freight 


3/ 

Third  Auxiliary  Model— 

Unmodified 

Modified  for 

Locational 

Interference 

20,6 

25.3 

5 

4 

53 

55 

42 

41 

100 

100 

Table  2 

Kesults  From  the  Third  Auxiliary  Model 
(data  for  consumer  products  manufacturer) 


gives  a tolerable  prediction  as  to  the  particular  warehouses  which  should 
be  selected  for  use:  of  the  top  ranking  12  (24)  warehouses  according  to 

the  model,  11  (14)  appeared  in  the  optimal  solution  of  the  integer  program) 

3/  Actually,  a slightly  different  version  of  the  outbound  transportation 
~ cost  term  of  (5)  was  used.  Instead  of  using  average  values  for  tj^  , 

which  were  not  conveniently  available,  a regression  was  performed  to 

- i 2 

R "*  ¥ 2 

fit  an  e;^ression  of  the  form  aD,  in  place  of  . 377  t.  p D 

K K K 

(note  that  a is  teiken  to  be  independent  of  k). 
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Finally,  consider  again  the  multiplicative  factoring  by  j of  all  in- 
bound transportation  rates,  the  change  which  foiled  the  second  auxiliary 
model  by  causing  an  increase  of  4 in  the  number  of  warehouses  open  in  the 
optimal  integer  programming  solution.  The  third  auxiliary  model  yields 
the  surprisingly  accurate  prediction  that  the  increase  will  be  4.4. 

Since  we  are  really  seeking  insights  rather  than  numbers,  it  is  even  more 
satisfying  that  an  intuitively  satisfactory  explanation  can  be  given  for 
the  increase.  The  explanation  is  straightforward  but  would  require  going 
into  more  numerical  details  than  seems  appropriate  here. 


- 21  - 


III.  CONCLUSION 


We  have  chosen  to  illustrate  the  general  methodological  approach  ad- 
vocated in  Sec.  I by  giving  more  than  one  auxiliary  model  for  the  same 
application.  This  was  done  to  convey  the  potential  diversity  of  auxiliary 
models  likely  to  be  available,  and  to  stress  that  developing  auxiliary 
models  can  be  an  iterative  design  process  guided  by  the  kind*  of  insights 
desired.  Such  models  cem  sometimes  be  taken  from  the  existing  literature 
on  analytical  solutions  to  highly  simplified  cases.  Sometimes  available 
amalytic  models  must  be  modified  or  new  ones  created  even  to  the  point  of 
losing  a closed  form  solution.  Of  course,  one  should  stop  short  of  the 
point  where  the  loss  of  auxiliary  model  tractability  due  to  added  com- 
plexity destroys  its  comparative  advantage  as  a source  of  insights. 

A rich  source  of  auxiliary  models  is  available  via  the  notion  of 
Lagrangean  relaxation  [2] . This  technical  device  is  aimed  at  producing 
highly  tractable  simplifications  of  difficult  mathematical  progrcimming 
problems.  Although  originally  conceived  for  purely  algorithmic  uses,  it 
turns  out  that  Lagrangean  relaxations  are  frequently  amenable  to  natural 
managerial  interpretation  and  hence  are  attractive  candidates  for  adop- 
tion as  auxiliary  models.  Initial  computational  experience  along  these 
lines  indicates  that  the  idea  will  be  a fruitful  one.  A very  recent 
paper  [4]  takes  the  concept  a step  further  auid  proposes  that  auxiliaty 
models  based  on  Lagrangean  relaxation  be  used  to  help  guide  the  aggre- 
gation choices  for  the  full  scale  model. 

It  seems  not  at  all  unlikely  to  me  that  auxiliary  models,  some  based 
on  analytic  prototypes  and  others  on  Lagrangean  releixation  or  other  form- 
alisms, will  eventually  be  common  in  strategic  applications  of  mathe- 
matical programming.  They  can  be  valuable  aids  for  designing  full  scale 


REFERENCES 


liJ 


1.  Bos,  H.C. , Spatial  Dispersion  of  Economic  Activity,  Rotterdam 
University  Press,  1965. 


2.  Geoffrion,  A.M.,  "Lagrcingean  Relaxation  for  Integer  Programming," 
Mathematical  Programming  Study  2,  1974. 


3.  Geoffrion,  A.M.,  "A  Guide  to  Con^uter-Assisted  Methods  for 

Distribution  Systems  Planning,"  Sloan  Management  Review,  Winter,  1975. 


4.  Geoffrion,  A.M.,  "Using  an  Auxiliary  Model  to  Guide  the  Design  of  a 
Very  Large  Logistic  Planning  Model,"  Discussion  Paper  No.  60,  Graduate 
School  of  Management,  UCLA,  April  1976. 


5.  Ignall,  E.J.,  P.  Kolesar  and  W.E.  Walker,  "Using  Simulation  to  Develop 
and  Validate  Analytical  Emergency  Service  Deployment  Models,"  P-5463, 
Rand  Corporation,  Santa  Monica,  August  1975. 


6.  Learner,  E.E.,  "Locational  Equilibria,"  Journal  of  Regional  Science, 
8,  2 (1968)  . 


7.  Magee,  J.F.  , Industrial  Logistics,  McGraw-Hill,  1968. 


8.  Meisel,  W.S.  and  D.C.  Collins,  "Repro-Modeling:  An  Approach  to 

Efficient  Model  Utilization  euid  Interpretation,"  IEEE  Trans,  on 
Systems,  Man,  cuid  Cybernetics,  Vol.  SMC-3,  No.  4,  July  1973. 


■» ‘ii  • -* — - ~-- 


I 


4 


